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Abstract 

. In this letter, we study the late-time evolution of a torsion cosmology only with 

the spin-0 + mode. We find three kinds of analytical solutions with a constant affine 
q ' scalar curvature. In the first case, it is not physical because the matter density will 

be negative. In the second case, it shows that the dark energy can be mimicked 
in the torsion cosmological model. In the third case, the characteristic of late-time 
evolution is similar to that of the universe of matter dominant. And we also find 
a kind of expression with the non-constant curvature that the periodic character 
of numerical calculation is only the reflection of solution in a specific period of 
f — \ evolution. Using these expressions, we shall be able to predict the evolution over 

the late-time. From this prediction, we know the fate of universe that the universe 
would expand forever, slowly asymtotically to a halt. 
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1 Introduction 



The cosmological revolutions of the past quarter century have changed ev- 
erything about our understanding of the future of universe. The current ob- 
servations, such as SNela (Supernovae type la), CMB (Cosmic Microwave 
Background) and large scale structure, converge on the fact that a spatially 
homogeneous and gravitationally repulsive energy component, referred to as 
dark energy, accounts for about 70 % of the energy density of universe. Some 
heuristic models that roughly describe the observable consequences of dark 
energy were proposed in recent years Dark energy can even behave as a 
phantom which effectively violates the weak energy condition [2]. When the 
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parameter of equation of state w is a constant less than —1, the universe ends 
up with a big rip singularity which is characterized by the divergence of curva- 
ture of the universe after a finite interval of time [3] . For the quintessence jl] 
and phantom [5] models, the future course of evolution is shown to critically 
depend on the potential. Specific phantom field models may be proposed to 
avoid the cosmic doomsday of big rip [6 J . 

On the other hand, it is inspiring that one can replace physical field by a 
geometry quantity in the dark energy model. The PGT (Poincare Gauge 
Theory of gravity) has a priori independent of local rotation and transla- 
tion potentials, which correspond to the metric-compatible connection 1-form 
pjw _ Y^dx a and orthonormal coframe ^ = e^dx a , where the metric is 
g = — <g> $° + Sij d 1 <g> $ and //, u, p • • • are Ad coordinate (holonomic) indices 
and i, j, k ■ ■ ■ are 3d. The gauge vector potential associates with field strength. 
In our case, they are curvature and torsion. PGT has been regarded as an in- 
teresting alternative of GR (general relativity) because of its gauge structure 
and geometric properties [7]. The bouncing cosmological model with torsion 
was suggested in Ref . [8] , but the torsion was imagined as playing role only at 
high densities in the early universe. Goenner et al. made a general survey of 
the torsion cosmology [9], in which the equations for all the PGT cases were 
discussed although they only solved in detail a few particular cases. Recently 
some authors have begun to investigate torsion as a possible reason of the 
accelerating universe [TO] . 

There are six possible dynamic connection modes [UJ, carrying certain spins 
and parity: 2 ± , 1 , . Some investigations showed that O 1 * 1 may well be the 
only acceptable dynamic PGT torsion modes [12] . The pseudoscalar mode 0~ 
is naturally driven by the intrinsic spin of elementary fermions, therefore it 
naturally interacts with such sources. Consequently, it is generally thought 
that axial torsion must be small and have small effects at the late time of 
cosmological evolution. This is a major reason why one does not focus on 
this mode at the late time. On the other hand, the scalar mode + does not 
interact in any direct obvious fashion with any known type of matter |13j . 
therefore one can imagine it as having significant magnitude and yet not be- 
ing conspicuously noticed. Furthermore, there is a critical non-zero value for 
the affine scalar curvature since + mode can interact indirectly through the 
non-linear equation. Nester and collaborators [14J consider an accounting for 
the accelerated universe in term of the dynamic scalar torsion. In Ref. [15] , it 
was shown that the dynamic Riemann-Cartan geometry could contribute an 
oscillating aspect to the acceleration expansion rate of the Universe. We have 
appeared analyzing the dynamics of this model [TB] . Applying the statefinder 
diagnostic to the torsion cosmology, we find that there are some typical char- 
acteristics [IT]. An extension model with the spin-0 + and spin-0~ modes was 
also considered, but the acceleration mechanism is still due to the spin-0 + 
mode [TS] . 
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In this paper, we study the late-time evolution of a torsion cosmology only 
with spin-0 + mode. In a constant affine scalar curvature case, we find three 
analytical solutions: the first is not physical, conflicting with the assumption 
of energy positivity; the second shows that the dark energy can be mimicked 
in the torsion cosmological model; the third shows that the behavior of late- 
time evolution is analogous to that of the universe of matter dominant. To 
satisfy the energy positivity requirement [12J, there is only a critical point 
(0, 0, 0) for the nonlinear system which is an asymptotically stable focus in 
the phase space (H, $, R) [16J. In non-constant curvature case, we also find a 
kind of expression showing that the periodic character of numerical calculation 
is only the reflection of solution in a specific period of evolution. Using these 
expressions, we shall be able to predict the evolution over the late-time. From 
the prediction, we know the fate of universe that the universe would expand 
forever, slowly asymtotically to a halt. 



2 The model 



Nester and collaborators [2] consider an accounting for the accelerated uni- 
verse in terms of a Riemann-Cartan geometry: dynamic scalar torsion. The 
torsion and curvature 2-forms are defined by 



r"E-ry a A/ = # + r^Ar (i) 

= ~RZp$ a A $ = dT» v + T£ A T pv (2) 
which satisfy the Bianchi identities, respectively, 



DT M = R p A 0", Di?^ = (3) 

Theoretical analysis of PGT led us to consider tendentiously two spin-0 modes. 
In this case, the gravitational Lagrangian density is 

1 3 («) h + \r 

C[»,r\ = -[- ao R+^a n T 2 + -R* + -E 2 ] (4) 

(n) 

where T is the algebraically irreducible parts of the torsion, R is the scalar 
curvature and E is the pseudoscalar curvature [18]. Note that a and d n are 
dimensionless parameters, and have the same dimension with i? -1 . 

Since current observations favor a flat universe, we will work in the spatially 
flat Robertson- Walker cosmological model. The isotropic orthonormal coframe 
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has the form: 



= t, #* = a{t)f (5) 
Because of isotropy, the only non-vanishing torsion tensor components are 

T; o = T] k = -2 X (t)e) k (6) 

where = e^k is the usual asymmetric tensor. Using the equation obtained 
by the variation with respect to the connection, we know the 0~ part couples 
to the axial spin vector of spin-- fermions, but + mode does not couple 
to any known source [TS]. Therefore, one can consider only spin-0 + mode. In 
other words, we only discuss the case of x{t) = 0. Furthermore, from the field 
equation one can finally give the necessary equations for the matter-dominated 
era to integrate (for a detailed discussion, see Ref.fIB] and we have made the 
replacement A Q — > — ao and A n — > 2a n , which is consistent with an earlier 
work [T6].) 



£ = -| + (9) 



H=J^R-1P^-2H 2 (7) 
6d2 6a2 

$ = ^ R _ l£n _ 3i y$ + I$2 (8) 
2a 2 2a 2 3 

where b = b + , \i = a 2 — ao, H = a/ a is Hubble parameter, and the energy 
density of matter component is 

K Pm = 1(J2 + ^)(3H - $) 2 - Ai? 2 - 3a 2 H* (10) 



The Newtonian limit requires a = — 1. 



3 The solutions of constant scalar curvature 



From Eq. (9), it is easy to find the scalar affine curvature remains a constant 
R = — 6/j,/b forever as long as its initial data has this special value [2]. In this 
case, Eq. (7) can be rewritten as 

4a 2 b 2 y J 
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The positivity of the kinetic energy requires a 2 > and b > [14], so we have 
the solution 

3C H 
H(t) = Ctan[^(to -t) + arctan(— -)] (12) 



where ( = \ij 1 \j2a 2 b and Hq = if (to)- However, such a choice conflicts with 
the assumption of energy positivity in the R = —6fi/b case. 

If we audaciously relax the parameter requirement for positive kinetic energy, 
i.e., a 2 < —1 and fi < 0, this phantom scenario will turn out to be interesting. 
Now we have the solution 

H(t) = ^ + ffo)(e-#o)- 1 exp[3e(t-t )]-l 

U (e + J ffo)(e-^o)- 1 exp[3e(t-to)] + l 1 ) 



where £ = \ij ^—2a 2 b. When t tends to infinity, H{t) —> £, so that the dark 
energy can be mimicked in the torsion cosmological model with a constant 
affine scalar curvature. Using the dynamical analysis, we have pointed out 
that there is a late-time de Sitter attractor [16]. Note that the solution (13) is 
just corresponding to the de Sitter attractor. 

Especially, as a 2 = — 1, we have a solution 



H=^^- (14) 

^_ 3ff [2+(2 + 2ff t)^] 

(2 + 3^)[l + (2 + 3F t)^] 
R = (16) 

where 



2 2 / 3 (-3#o + $ c 
3#o - 2$ n 



(17) 



From Eq. (14), we have 

.2 + 3H t 2 
a = ao{- — — — — J 3 
2 + 3Hoto 



Note that the behavior of the late-time evolution is analogous to that of the 
universe of matter dominant. In Fig. 1, we plot the rajectories in the phase 
space H-$ when we take a 2 < —1, b > and R = —Qji/b. Obviously, if a 2 = 
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Fig. 1. The rajectories in the phase space H-$> fixing 02 < —1 in 1(a) and 02 = — 1 
in 1(b), respectively, in R = —6(j,/b case. 

or —1, the slightest change in 02 leads to a radical change in the behavior of 
solutions. Therefore, we have a bifurcation at a 2 = or — 1 in Eqs. (7)- (8) and 
R = -6/V&. 



4 The solution of non-constant scalar curvature 



The numerical analyses show that H, <£> and R have a periodic character at 
late-time of the evolution for 02 > and b > 0, approximatively [15]. Using 
the dynamical analysis [16J and statefmder diagnostic [17], we find that this 
character is corresponding to an asymptotically stable focus. In Fig. 2, we 
plot evolving trajectory in H-Q-R space, where we have chosen a<i = 2 and 
b = 2/£q- We find easily that the evolving trajectories tend to the focus (0,0,0). 
If we consider the linearized equations, then Eqs.(7)-(9) can be reduced to 



H=-^R, i> = ^-R, R = ~^ (19) 
oa 2 2a 2 

The linearized system (fT9l) has an exact periodic solution 
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Fig. 2. We illustrate evolution trajectory in the phase space, where have chosen 
02 = 2, b = 2/tg. The surfaces are the nullclines of nonlinear dynamical system, and 
the curve is the evolution trajectory. It shows clearly that at late time the phase 
line has a oscillation feature. 

H=-aR sinut + cos cut + H - 
$ = — /3 _1 i? sin ut + $ cos ut 

R = R cos ut + /3$o sin ut (20) 

where u = y^g, a = y^ , /3 = ^ and # = #(0), $ = $(0) and 
i?o = -R(O) are initial values. Obviously, (if, 0,0) is a critical line of center for 
the linearized solution. In other words, there are only exact periodic solutions 
for the linearized system, but there are quasi-periodic solutions near the fo- 
cus for the coupled nonlinear equations. This property of quasi-periodic also 
appears in the statfmder diagnostic with the case of ci2 > [TTj . According to 
nonlinear equations (7)-(9), we can obtain the critical points and study the 
stability of these points. To study the stability of the critical point (0,0,0), 
we write the variables near (0, 0, 0) in the form H = AH, $ = A$ and 
R = AR, where AH, A<3> and AR are the perturbations of the variables near 
the critical point (0,0,0). Substituting the expressions into Eqs. (7)-(9), we 
obtain the corresponding eigenvalues (0, —y—^i \[H^b)- Therefore, the criti- 
cal point (0, 0, 0) is a stable focus when a 2 > and b > 0. We illustrate the 
evolution trajectory in the phase space (H, <3>, R) in Fig. 2, where have chosen 
a,2 = 2, b = 2/io- Nullcline surfaces in Fig. 2 divided the phase-space (H, $, R) 
into some domains, in each of which there is definite behavior of the evolution 
trajectory. It is obvious that the evolution trajectory helically tends to the 
critical point (0, 0, 0) after crossing the nullcline surface. 

For the late-time behavior of non-linear system Eqs.(7)-(9), we should de- 
termine the behavior at infinity. We substitute r = \ in Eqs. (7)- (9), and 
show the Laurent expansions around the point r = 0. Clearly, as t — > oo, 
r — > 0. Thus, we are interested in the behavior of equations for r = 0. The 
non-linear system clearly shows that the solutions have regular behavior at 
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t = 0, and have the Taylor expansions: H(t) = E£U h n r n , $(t) = Y^=?, < PnT n 
and R(t) = Y^=i r rJ n ■ So we find an approximate formula up to t~ N order as 
follows 



N h 



n=l 
N 



t N 



\J3n sin cut + 7tv cos cut] 



1 ,3(5 



A? 



37. 



A? 



n=3 
N 



R(t) = E 



n=2 



1 r -127iv 
t N[ uob 



cos cut] 



sin cut H — cos cut 

cuo 



(21) 
(22) 
(23) 



where cu = \J~^i and h n , (p n and r n are undetermined coefficients. The coeffi- 
cients /3jv and 7a? can be written as the terms of initial values. 

From a theoretical point of view, the Eqs. (21)-(23) are equivalent to the 
Laurent expansions around the point r = when iV — > oo. From a practical 
point of view, the Eqs. (21)-(23) are more effective than the expansions of 
minus-power terms because of the latter only determined by the coefficients 
hi, ■ • ■ , h at— i j tpi,---, (pN-i\ r ii ' ' ' ; fN-i for a definite N. Therefore, the former 
is an approximate formula up to t~ N order and the latter is one up to t - ^ -1 ) 
order. For example, we take N = 3 and the approximate solution up to t~ 3 
term is 



rr (f) _ 2 1 ^0 U 

tfi u A 4 ^ n A 4 n °°b ub . . 



+(-^B„S„ + ^S + ^4> C - ^-C„) cos^f] (24) 
12 ytg o ytg 



+ (-^i?o^o + 7^0 + $ oC - ^G)) coscut] (25) 

4/X 3/xtg 3tg 



^ " 3^ + ^ [ lo" " 3tjj " 3^ ° + 9^ 

+ Wo - xUo - ^ oC + ^rC ) sincut 
3tg cub 3cobt 



+ (^ S - -^-S + R Co - ^2 Cg) cos cut] (26) 
cob 3cobt 3tg 



where H = if(t ),<3> = $(to) and i? = R(t ), and So = sincuto,Co = 
coscut . Obviously, H{t) ->■ 0,$(t) ->■ and i?(t) ->■ 0, when t tends to 
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3(b) 



Fig. 3. In the non-constant curvature case, we plot the behavior of late-time evolu- 



tion. In 3(a), we have chosen a 2 = 1,6 



and to = 20; In 3(b), we have fixed 



a 2 = 1, 6 



8000 

TT 2 t 



¥ and t = 20. 







infinity. Therefore, the fate is that the universe would expand forever, slowly 
asymtotically to a halt. 

Next, we give the reason why some numerical calculations look as periodicity 
as mentioned in Ref . [15] . To put it bluntly, this is a reflection of the late-time 
solution at specified period of evolution. If we take t = to + At, At G [0, 5} and 
^ <C 5 <C toi Eqs. (24)-(26) can be reduced to a solution, whose behavior looks 
like a periodic solution of (20) of the linearized system. Taking the parameter 
values, we plot Hubble parameter H(t) at the late-times in Fig. 3. 

If we take the approximate expression of N = 2, when t > 10, 10 2 and 2 x 10 2 , 
the error A = ferggg=al < 0.O66, 0.012 and 0.002, respectively. Therefore, 
if we want to obtain an effective approximate expression, N is required to 
be large enough. In Fig. 4, a direct quantitative connection is made between 
analytical expressions and the numerical solutions presented by Shie et al. 
H5]. We take initial conditions H(l) = 1,$(1) = 1.4,^(1) = 1.53, /1 = 1.09 
and b — 1.4 which are the same as those in Ref. [T5] . 
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Fig. 4. In the non-constant curvature case, we plot the behavior of late-time evolu- 
tion H(t) for our analytical expression and numerical solution obtained by Shie et 
al., respectively. Black line represents the numerical solution (t < 3), and red line 
represents the analytical solution (JV = 50). 

5 Conclusion 



In this paper, we study the evolution of a torsion cosmology only with the 
"scalar torsion" mode + . This mode has some distinctive and interesting 
qualities. It can be considered as a "phantom field", because it does not in- 
teract with any known matters directly. It only interacts via the gravitational 
equations. Therefore, + mode can drive the universe to accelerate at present. 
However, we have to investigate the late-time solution if we want to know 
the fate of universe. We find a kind of late-time solution which tends to the 
focus (0,0,0) in the phase space(if, $, R). Under the solution , the feature of 
evolution is universal with generic choices of the parameters. Distinctive and 
interesting results are: 

i) In the constant affine curvature case, there are a kind of non-physical 
solution and two kinds of physical solutions at late-time. One of physical 
solutions corresponds to the de Sitter attractor, and the other is analo- 
gous to the universe of matter dominant. In the view-point of dynamical 
analysis, 02 = and —1 are the bifurcations. 

ii) In the non-constant affine curvature case, we find a kind of expression 
corresponding to the universe which would expand forever, slowly asymp- 
totically to a halt. The specific behaviors of late-time evolution differ from 
those at different periods. If we considered the specified period of evolu- 
tion t = to + At, At G [0,5], and satisfied — <ti 5 <ti to, the behavior 
looks like a periodic solution. 

Furthermore, if we want to establish a realistic cosmological model, we have to 
carry out the solutions of t < t period, and compare these analytical solutions 
with observations to determine the model parameters. For example, using SN 
la data, we constrain the parameters of torsion cosmology and get their best 
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values. Then we can discuss the properties of the realistic cosmological models. 
These issues will be considered elsewhere. 

Recently, the cosmological model with even and odd parity modes in PG has 
been considered by Baekler, Hehl and Nester [19J. They extended the parity 
violating quadratic Lagrangian VL to the general gravitational Lagrangian 
V± = V- + V+, where V + is the parity conserving quadratic Lagrangian. This 
model generalized the torsion cosmologies which were presented by Nester et 
al. [T2lll4lll5|IT8] . Next, we will set about investigations for dynamics of this 
cosmological model. 
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